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Abstract 

Various QCD correlators are calculated in the instanton liquid model in zero- 
mode approximation and 1/Nc expansion. Previous works are extended by including 
dynamical quark loops. In contrast to the original "perturbative" 1/Nc expansion 
not all quark loops are suppressed. In the flavor singlet meson correlators a chain of 
quark bubbles survives the A'^c ^ oo limit causing a massive r]' in the pseudoscalar 
correlator while keeping massless pions in the triplet correlator. The correlators are 
plotted and meson masses and couplings are obtained from a spectral fit. They are 
compared to the values obtained from numerical studies of the instanton liquid and 
to experimental results. 
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1 Introduction 



The fundamental theories which describe the observed interactions are all gauge theories. 
They describe gravitation, electroweak and strong interaction. In the quantized version 
the forces between particles (fermions) are mediated by gauge bosons. QCD describes 
the strong interaction between quarks and gluons. It is a SU(3) gauge theory with the 
LagrangianQ 



Taking into account only the light quarks u, d and sometimes s and setting their mass 
to zero one arrives at a theory with only one parameter, the gauge coupling constant g, 
which is actually no parameter because of dimensional transmutation. Despite the fact 
that QCD looks such simple (e.g. compared to electroweak interaction) it is very hard to 
solve this theory due to nonperturbative effects. To clarify the role of instantons in QCD 
let me first give a list of the most important methods used to tackle QCD starting with 
very general methods applicable to any quantum field theory (QFT) and ending with 
more specific approaches: 

• Axiomatic field theory: Wightman/Oswalder&Schrader have stated a set of 
Minkowskian/Euclidian axioms for vacuum correlators a general QFT should re- 
spect (analyticity, regularity, Lorentz invariance, locality, ...). It is clear that the 
theorems derived from these axioms have to be very general because no Lagrangian 
is used 0. 

^ Throughout the whole work the euchdian formulation of QCD is used pl| . 
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Haag-Ruelle/LSZ Theory: S-Matrix elements are related to vacuum correlators. 
The S-Matrix is the central object for particle phenomenology containing such im- 
portant informations as cross sections, form factors, structure functions, .... For 
theoretical studies vacuum correlators are more suitable because they avoid the need 
of explicitly constructing the Hilbert space which is an extremely complicated task 
beyond perturbation theory. 

Quantization: One might think that quantization should not appear in a list of 
methods for solving QFT because it is the method to obtain and define a QFT. On 
the other hand, besides canonical quantization there are other ways of quantizing a 
theory. The most popular is the path integral quantization . Usually in textbooks 
for particle physics it is only used as an abbreviation to derive theorems more quickly. 
In general (beyond perturbation theory), different quantization methods lead to 
different physical and mathematical insights and different methods to solve the 
theory. Variants of the path integral quantization are the random walk quantization 
used in lattice theories and stochastic quantization. 

(Broken) Symmetry: Every degree of freedom like spin, flavor and color is the 
possible origin of (approximate) symmetries like SU{2Nf) or subgroups and SU{3) 
gauge invariance. Conserved currents and Ward identities can be obtained. In the 
case of light quarks, one further has approximate chiral symmetry leading to PCAC, 
axial Ward identities, current algebra theorems, soft pion physics,. . . . Furthermore, 
QCD with massless quarks possesses an anomalously broken scale invariance which 
is the origin of the huge field of renormalization group techniques 0. 

Perturbation theory: Due to asymptotic freedom the coupling constant g decreases 
for high energies and perturbation theory in g is applicable. QCD can thus be 
solved for processes which involve only momenta of say more than 1 GeV. The 
small distance behaviour of vacuum correlators is thus calculable {x < 0.2 fm). 

Operator Product Expansion (OPE): An improvement of perturbation theory is 
to separate the small distance physics from large distance effects. The former is 
contained in the so-called Wilson coefficients calculated perturbatively. The latter 
nonperturbative effects are contained in a few vacuum or hadron expectation values 
of local operators which have to be determined from phenomenology or uncertain 



assumptions [|TT|. Vacuum correlators can be obtained up to distances of x < 
0.3. . .0.5 fm. 

QCD Sum rules: QCD sum rules are widely used to determine hadron masses and 
couplings. The general method is to assume the existence of certain hadrons and 
take a resonance+continuum ansatz in the Minkowskian region for some correlator. 
The Euclidian correlator is calculated from theory (OPE, lattice, instantons). Via 
dispersion relations one may match both in some Euclidian window by fitting the 
hadron parameters which leads to a prediction for them. 

Effective Theories: One may construct effective Lagrangians containing mesons 
and/or baryons more or less motivated by QCD or history or other physical branches. 
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Their parameters are determined from experiment or QCD as far as possible. 



A variety of phenomena have been explained qualitatively and calculated quantitatively 
with the methods listed above, but the large distance behaviour of QCD is still unsolved. 
There are at least two problems belonging to this domain: chiral symmetry breaking and 
confinement. Up to now chiral symmetry breaking was assumed and the consequences 
such as Golstone bosons were discussed within this assumption. In OPE one takes the 
nonzero values of the quark and other condensates from experiment but has no possibility 
to predict them from theory. The quark condensate is the order parameter of chiral 
symmetry and a nonzero values indicates spontaneous breaking of this symmetry (SBCS). 
Confinement has also to be assumed. 

These two problems are not so unsolved as it might appear. Up to now I have not 
mentioned two further approaches: 



Lattice QCD: In principle the method is very simple. The continuum is replaced by 
a fine lattice covering a large but finite volume. The path integral is thus replaced 
by a finite number of integrals evaluated numerically. All vacuum correlators can 
be obtained for arbitrary Euclidian distances. Confinement and SBCS have been 
shown and other hadronic parameters are obtained. In practice, lattice calculations 
are much less straightforward than this scetchy description might suggest 0. 

Instantons: As in lattice QCD one evaluates the Euclidian path integrals but 
now in semiclassical approximation. In addition to the global minimum of the QCD 
action = used by perturbation theory there are many other local minima called 
instantons which have to be taken into account. Inclusion of light quarks leads to 
an effective 2Nf quark vertex responsible for SBCS. Although confinement cannot 



be explained, a lot of hadron parameters can nevertheless be calculated pO| , ^ 
suggesting that confinement is not essential for the properties of hadrons. 



In this article I want to give a quantitative and systematic study of the implications of 
instantons to hadron properties. 

Chapter |^ is an introduction to the semiclassical evaluation of nontrivial integrals. After 
developing the method in the finite dimensional case the partition function of QCD is 
considered and the instanton liquid model is introduced. 

In chapter |] the propagator of a light quark is calculated. The approximations which 
have to be made are stated and discussed. It is shown that for one quark flavor the 
1/N(. expansion is exact. The constituent quark masses and the quark condensates are 
calculated for u,d and s quarks. 

The same approximations are used in chapter ^ to calculate the 4 point functions. Special 
attention is payed to the singlet correlator where a chain of quark loops contributes and 
is not suppressed in the \aige-Nc limit. Within one and the same approximation we get 
Goldstone bosons in the pseudoscalar triplet correlator but no massless singlet boson. 
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In chapter || meson correlators are discussed and plotted. Employing a spectral ansatz 
it is possible to extract various meson masses and couplings. They are compared to the 
values obtained from extensive numerical studies of the instanton liquid [21| and to the 
experimental values. 



Conclusions are given in chapter 



2 Theory of Instanton Liquid 

Separating Gaussian from non-Gaussian degrees of freedom 

In the path integral formulation solving QCD or any other QFT is equivalent to calculating 
the partition function Z including external source terms. Of course this is not a simple 
task because one has an infinite number of degrees of freedom. One way to tackle this 
problem is to reduce the number of degrees of freedom by integrating the (nearly) gaussian 
degrees perturbatively leaving the non-Gaussian degrees, called collective coordinates, to 
be treated with other methods. 

Let me first describe this method most generally without refering to QFT. Our aim is to 
calculate 

Z= tee-^["l (2) 

where S is some real valued function depending on the n dimensional vector x G iR". Let 
the integral be dominated by values of x which lie in the vicinity of a k dimensional sub- 
manifold of iR" which may be parametrisized by x = /(7), 7 G iR^. Vectors in the tangent 
space of /(7) are called (approximate) zeromodes, because 5'[a;(7)] is (nearly) constant in 
this direction. Usually /(7) represents some degenerate or approximate minima of 5*. In 
addition /(7) may contain points which do not contribute much to the functional integral, 
this will cause no error, but /(7) must not forget any significant points. In figure a two 
dimensional example is shown containing a river (one dimensional manifold) with steep 
mountains aside his banks. So the integral will be dominated by the shaded area. This 
area can be parametrized in the following way: 

x = f{j) + y , x,yelR^ , 7 G iR'^ . (3) 

To make this representation unique, we have to demand k linear extra conditions to y: 

y-9iil) = Ci , l<i<k (4) 

E.g., Cj = and (71(7) = df{^)/d'^i would fix y to be orthogonal to the "river" or equally 
stated: It disallows fluctuations in the non-Gaussian zeromode direction. Now every 
point X (at least in the vicinity of the river) can uniquely be described by y satisfying the 
condition @) and by 7 via (^. All we have to do now is to represent Z in terms of the 
fluctuation vector y and the collective coordinates 7. A convenient way is to introduce a 
Faddeev-Popov unit 



1 = |rf^7 d^y5\yg,{^) - c,) 5"(/(7) + y - x) $(x) 



(5) 
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which serves as a definition of $. Inserting this into and integrating over x yields 

Z = |rf^7A5'(2/-^7.(7)-Q)<f(/(7) + y)e-^[^(^)+^l . (6) 
The ^/-integration in can trivially be performed: 

<|.-i(x)= /dV5^(x-/(7')k.(7')-Q) . (7) 



For X = /(7) + y the 5-function in the integral only contributes for 7' = 7. Thus we can 
expand the (5-argument up to linear order in 7' — 7 

= \detij{y-djgi- gi-djf)\~'^ ^ 
with y-gi{-f) = q . 

Inserting $ into (H) we get 

Z = J d^^<ry5\yg,{^) - q) det,, {yd.g^i) - gAl)-Qof^l)) e-^'^^^^+^l • (9) 

One may write Z in a slightly different form, usually used in QFT, because it is more 
suitable for semiclassical approximations. Z is independent of q and therefore, although 
the r.h.s. of (P) explicitly contains q, it is actually independent of it. We can smooth the 
5-function by a further multiplication with 

1 = (27rO"^/'yc/^ce-i^'"' . (10) 

The determinant can be written in the form 

detA= fdr]dr]e'^^'' (11) 



where rj are anticommuting grassmann variables (ghosts). Inserting (^) and ([11]) into (^) 
and performing the c-integration we finally get 



s,f = Y^y'^(J:9^{l)9nlyjy , (12) 

Sfpg = J2^i(y-^j9iil) - 9i{l)-djf{l))r]j . 



^In the following we will omit the absolute bars; when necessary they have to be reinstated. 
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Effective QCD Lagrangian in a background field 

Now it is time to return to QCD 

Z = jDA.e-^-^'^^^ , SyM[A]= Jdx^^G^^G^'' . (13) 

The background configurations whicli (approximately) minimize Sym will be denoted by 
^^(7). A general gauge field can be written in the form 

A, = {A,{^) + B,f , 7eiR' (14) 

where 5^ are fluctuations around the background and A^ is a gauge transformed fleld 

A^ = SA^S^ + tSd^S^ , S = e'^ e SU{N^) . (15) 

As in the flnite dimensional case we have to make this representation unique by introducing 
extra conditions, 

D^{A)B^{x)=Cix) , (16) 
to flx the gauge of the fluctuating fleld and 

Jd^x^l{x;^)B^{x)=Ci , l<i<k (17) 

to avoid fluctuations in (approximately) zero mode direction. The derivation of an effec- 
tive action similar to (|T2]) can now be performed in full analogy to the previous case with 
only some notational complication. There is the following correspondence: 



flnite example : i x y "ji gi 
QCD : i,x A B -fi,n{x) 

The Faddeev-Popov unit has the form 



d>'^DnDB^5iD^iA)B^)5''i /#S^rf^x)5((4 + B^f - A^)^A^] 



(19) 



The steps to get an expression for $ are now: 

Add primes to 7, f2 and B, insert A = A + B, linearize the last 5- argument around 
B' = Bai perform the functional B' integration and linearize the remaining 5 arguments 



around 7' = 7 and Vt' = 0. Omitting the details of this calculation one gets [jl9 



$-^(A(7) + B) = j ct'i Dn'5^{Xi) 5{Y) 



(20) 



dA^ d^^jj) 



BA{i^-l,)+i,'^D,{A + B)n' 



(21) 
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- 8A 

y = T.DM + B)^{i-l,)+D,{A)DM + B)^' . (22) 
From ([16|) , (|17D , (|2ll) and (^) one can read off the form of tlie partition function Z: 

Z = N{^) Jd'''yDB^Dr]Df]6\J'^'^B^d^)e-^Q''''^^'^''^'^^ , (23) 

Sqcd = SyM[A + B]~Sgf[A,B] + SFPG[A,B,7],r]] , 
Sym = jd'x^^Gl^G>:^{A + B) , 

: ^Jd'x{D,{A)B,r , (24) 



Sfpg = 



d^^f^iV^ rf 



1j 

+ Yl Jd'xn^^;D^{A + B)r]{x) + Jd'xY,nix)D^{A + B) 

i j 

+ Jd^xf]{x)D^{A)D^{A + B)ri{x) . 



Sqcd does not depend on tlie gauge parameter Q. For this reason the Q integration can 
be absorbed in the normalization factor N{^). r]{x) are the usual ghost fields originating 
from the gauge fixing. For every extra condition (|17D one gets an additional ghost variable 
r]i. For A = and no extra condition [k = 0) the action given above just reduces to the 
usual QCD action including Faddeev-Popov ghosts in gauge 

Sgf = ^Jd^x{d^B^y ' Sfpg = Jd^xr]{x)d^D^{B)r]{x) . (25) 

Note that the action (^) is still exact with the non-harmonic degrees of freedom 7, now 
separated from the hopefully more gaussian ones, B^ and 77. 



For small coupling g it is now possible to establish feynman rules from (|24D in analogy 
to the case with no background. For this one has to know the "free" gluon, ghost and 
quark propagator in a given background A. If A is a non constant field even this is a 
very complicated task in contrast to usual perturbation theory around A = 0. For a 
multi-instanton configuration explicit expressions for the gluon and ghost propagator are 



derived in [|T6l , ^ and for light quark propagators in chapter | 
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The Semiclassical Limit 

Before developing perturbation theory to all orders it is wise to study the semiclassical 
limit where one keeps only terms up to quadratic order in the fields. In QCD (and many 
other field theories) this is equivalent to lowest order perturbation theory, but around a 
very nontrivial background! 

Up to now we have not specified ip^^. A natural choice would be ip^^ = dA^/d'^i to fix the 
fiuctuations to be orthogonal to the zero modes. Somewhat more convenient is to bring 
ip^^ in background gauge: 

(BA \^ 

= (^-^J with n such that D^{A)i)'^ = . (26) 

Furthermore we assume that A minimizes the gauge action Sym which is true for widely 
separated instantons thus neglecting linear terms Sqcd- 

Up to quadratic order in the fields one has 

Sqcd = Sym[A]+ Jd^^B^K^,{A)B, + Jd'xr]{x)D\A)r]{x) 

BA r - BA 

+ T.V^^;^m+ d'^T.vix)D,{A)-J^V,+0{^^^^^ , (27) 



1, 

Performing the integration over gauge fields and ghosts one gets an effective action de- 
pending only on the collective coordinates 7^: 



d^^e-'^^"^^^ (29) 

The 5-function in (|23|) causes a restriction of the gauge field fluctuation to be orthogonal 
to t/'^. K'^j^ is defined as K^^, projected to the space orthogonal to ipli, Det' takes into 
account all eigenvalues of K'^^^ except the k zeromodes caused by the projection. 

The Instanton Gas Approximation 

For well separated pseudoparticles A{x) = J2Ai{x) the partition function Z factorizes 
like 

^=^n^. (31) 
• 1=1 
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where Zj is the partition function in a one pseudoparticle background, i.e. with Aj inserted 
in (^) instead of A. It was a great deal to evaluate the functional determinants in the 
background of one instantonQ 



li 



{zi,Oj,Pj,Qj) 



oTvL 



2p2 



X Zi 



(32) 



(location, orientation, radius, topological charge) 



To get finite results Zi has to be normalized to the A" = case, regularized and renor- 



malized. The final result of the very complicated calculation is W\ |T3 

1 r poo 

- J2 Jd%dOidpiD{pj) = V,j^ dpD{p) 

Dip) = 



reg 



1 4.6e-i-6^9iv. 



5'o(p)'^^e-^«(^) 



p5 7r2(iV,-l)!(iV,-2)! 



(33) 



In^ 



A = A 



PV 



D{p) is the instanton density, g{p) the running coupling constant, b = and b' = y-^c- 



Quarks 

Additional fields coupled in a gauge invariant way to the gluon field can simply be incor- 
porated by adding the appropriate lagrangian with gauge field A replaced hj A + B and 
performing the functional integration over the new fields. So every quark contributes an 
extra factor 

'D^Dfre-/'^"*^'-^*™)* = Det(z^+ im) (34) 



I' 



to the partition function Z where 



(35) 



is the covariant derivative. In the semiclassical approximation can be set to zero. -D(p) 
has to be multiplied by the fermionic factor 



F{mp) 

and b and b' are now 



1.34mp(l + m p In(mp) + . . .) for mp <^ 1 
+ . . . for mp ^ 1 



b = —Nc - -Nf , b' = —Nj. NcNf + ^ . 

3 3-^ 3^ 3 -^2 AT, 



(36) 



(37) 



^ To simplifiy notations we will treat instantons and anti-instantons on the same footing. Both will 
be called instantons and are distinguished by their topological charge Q/ = ± if necessary. 
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The Instanton Liquid Model 

The probabihy of smaU size instantons is low because D{p) vanishes rapidly for small 
distances. On the other hand for large distances -D(p) blows up and soon gets large. This 
is the origin of the infrared problem which made a lot of people no longer believing in 
instanton physics. Those who were not deterred by that have thought of the following 
outcome |^ . For larger and larger distances, the vacuum gets more and more filled with 
instantons of increasing size. At some scale the instanton gas approximation breaks down 
and one has to consider the interaction between instantons which might be repulsive to 
stabilize the medium. The stabilization might occur at distances at which a semiclassical 
treatment is still possible and at densities at which the various instantons are still well 
separated objects - say - not much deformed through their interaction. So there is a 
narrow region of allowed values for the instanton radius. This picture of the vacuum is 
called the instanton liqiud model. The idea has been confirmed in the course of years by 
very different approaches: 



Hardcore assumption |]ll 
Variational Approach p9|| 
Numerical studies 
Phenomenological success El 



The picture has now become generally accepted at least by those who believe in instanton 
physics and it seems that the vacuum can be described by effectively independent instan- 
tons of size p = 600 MeV~ and mean distance Lq = 200 MeV. The integral instanton 



density is fixed by the experimentally known gluon condensate p5 



n = N/V, = l/Ll = ^ < >= (200 MeV)^,^. (38) 

The ratio Lq/ p is estimated in different works to be 

{Lolp)theor. = 3.0 . . . 3.2 (39) 



3 Light Quark Propagator 



In this chapter the average quark propagator in the multiinstanton background will be 
calculated. In the first section the instanton background is treated as a classical external 
perturbation, but the background field is not small (e.g. in the coupling ^f) and therefore 
we have to sum up all feynman graphs. This is possible in the case of one quark flavour 
within the so-called zeromode approximation. The quark condensate and a constituent 
quark mass are extracted from the quark propagator. In the last section it is shown that 
the case of two or more quark flavours can be reduced to the one flavour case in the limit 
of Nc — > oo. The results obtained in the one flavour case are therefore still valid when 
making this further approximation. 
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Perturbation Theory in the Multi-instanton Background 

It is well known how to calculate correlators in the presence of an external classical gauge 
field at least as perturbation series in powers of the external field A^(x). In the case of 
QCD (or more accurately in classical chromodynamics) within the instanton liquid model 
the external field is a sum of well separated scattcrcrs A = J2i Aj called instantons with 
a fixed radius p and distributed randomly and independently in Euclidian space. 

For a while we will restrict ourself to the case of one quark fiavor and ignore gluon loops. 
The Euclidian feynman rules have the following form: 

— So , 



1 



jjl + im 
Ai 

^ ^ 4i ■ (40) 



In operator notation, 



{p\So\q) = -t^Hp - q) , {x\4j\y) ^ Mx)5(x - y) , (41) 

~r tin 



graphs are simply alternating chains of Sq and Aj. To average a graph over the instanton 
parameters one has to perform the following integration for each instanton: 

i...h-^Jiv--^Z^^Jd%JdO,... . (42) 

For example 




= / d^id^j{p\So4iSo4iSo4jSo\q)Tr{So4jSo4i) ■ (43) 

The origin of the factor N'^ is the summation over all pairs of different instantons (/, J) 
yielding a factor N{N — 1) fa N'^. The quark loop is the origin of the minus sign and of 
the functional trace " Tr" . 
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Exact Scattering Amplitude in the one Instanton Background 

Perturbation theory is suitable to study scattering processes. To achieve chiral symmetry 
breaking or bound states one has to sum up infinite series of a subclass of graphs or solve 
Schwinger Dyson or Bethe Selpeter equations. 

The first thing we can do is to sum up successive scatterings at one instanton 

Ar X At )kA, 



^(©^ := + ^ ' + ■ ^ '1 ^ ■ +••• , 

Vi := 4i + 4jSo4i + 4iSo4iSo4i + ■■■ = So\Sj - ^o)5'o-' (44) 

where 

Sr = is,' - 4ir' (45) 
is the quark propagator in the one instanton background. 

Zero Mode Approximation 

The quark propagator in the one instanton background can be represented in terms of 
the eigenvalues Aj and eigenfunctions ipi of the Dirac operator 

(.^ - 4/)^. = A.V'. =^ {ASi\y) = Y. '^f^f-^^^ . (46) 

^ Xi + im 

There is one zero eigenvalue ii/'V/ = (tappendix ^ which makes the propagator singular 
in the chiral limit: 

S,JMM + Sf^^ . (47) 
im 

In the so called zero mode approximation one replaces the non zero mode part by the free 
propagator 



—45— =5.-5o^^^^ . (48) 

im 

Although this approximation is good for large as well as for small momenta it may be 
bad for intermediate ones, but what is more important is the fact that it is a wild ap- 
proximation and so might violate general theorems like Ward identities. In contrast, all 
other approximations we make are of systematic nature respecting all known symmetries 
of QCD. 

• semiclassical approximation (systematic) 

• multi-instanton background ("systematic") 

• large Nc expansion (systematic) (see next section) 

• zero mode approximation (wild) 

Note that every choice of a background gauge field is "systematic" in the sense of re- 
specting the symmetries of QCD as long as the background satisfies these symmetries on 
average. 

In the following sections we will see that the advantage of the zeromode approximation is 
so great that we cannot disregard this simplification. 
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Effective Vertex in the Multi-instanton Background 

Let us consider a quark line witli two scatterings at Vj and insert in between a number 
of instantons wfiicli differ from I and from all other instantons occuring elsewhere in the 
graph. This enables us to average over these enclosed instantons independently from the 
rest of the graph. Summation over all possible insertions with at least one instanton 
just yields the exact quark propagator minus the free propagator. Remember that direct 
repeated scattering at Aj is already included in Vj. 



Let us define 



^ := € ' = S-So 

(49) 

Mi = Vi + Vi{S-So)Vi + Vi{S-So)Vi{S-So)Vi + ... 
= Vi + Vi{S - So)Mi 

This equation can be solved for Mj with the following ansatz: 

Mi = -So'\4ji){HSo' (50) 
Inserting M/ and Vj into (^) we get 

^S^'\^|JJ){^J\S^' = —(1 + -{^|JI\S,\S - S,)S^'\^J))S,'\^|JJ){^J\S^' (51) 
=^ l^ = m + t{^j\So\S-So)S^'\^j) . (52) 



A Nice Cancellation 



It is possible to arrange the graphs in such a way that every Mj occurs only once. Consider 
a graph containing two scattering processes Mj at the same instanton. The interesting 
part of the graph has the following form 



P 



ap 


'^i{p)ip\{qy 


Olq Cer 




tp 




Iq If 





(53) 



ap/ip are the color/Dirac indices at the quark leg with momentum p. In a physical process 
in addition to the graph containing the above subgraph there exists another graph with 
only two quark lines interchanged. 

(54) 



ap 








tp 
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As usual, the interchange of two quark hnes causes a minus sign in the amphtude. In- 
specting the two expressions, we see that they coincide except for the sign, thus there 
exists a complete cancellation 




Whenever an Mj occurs twice or more than twice in a graph there exists another graph 
with opposite sign. Both contributions cancel each other and can be ignored. So a quark 
can scatter only once at every instanton. This can be seen in another way: Because of 
Fermi statistics every state can be occupied only once, and there is only one state for each 
quark in the zero mode approximation, namely the zeromode. 

There are two equivalent descriptions of feynman graphs: 

1. Draw all topologically distinct graphs with non-numerated vertices and assign a 
symmetry factor to each graph, 

2. Draw all topologically distinct graphs with numerated vertices and assign a factor 
1/V, where V is the number of vertices. 

If all possible graphs containing Mj are allowed it is not difficult to see within the second 
description that they can really be paired as stated above. 

Renormalization of the Instanton-Density 

Up to now the cancellation is incomplete because not all graphs are allowed. Consider 
e.g. 




not allowed not allowed 

As in the case of Vj both graphs are not allowed. Another example is 




+ —^^^^-^ = 

not allowed most general tadpole 

It is a general fact that all disallowed graphs can be paired with other disallowed graphs 
or with tadpoles and vice versa. 
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What we have to do is to "disallow" all tadpole graphs. Every Mj can be surrounded by 
tadpoles which contribute with a universal multiplicative factor which can be absorbed in 
a redefinition of the instanton density ur. Using this renormalized density nn the pairing 
is now perfect and the statement "every Mj occurs only once" becomes true. 

One further can show that in the presence of dynamical quarks this renormalized density 
has to be identified with the gluon condensate instead of the "bare" density because the 
same tadpoles contribute to the gluon condensate too. 



Selfconsistency Equation for the Quark Propagator 

Quark loops are no longer possible because they cannot be connected to another part of 
the graph via a common instanton. All graphs which can contribute to the propagator 
are chains of different MjS. 

— € — = ( — ' — + + -^^^y*^)-^ + — 

\ / i^j^... 



The M[s can be averaged independently 



M(p):=z(M,) = ^pV2(p) , (56) 

where ip' is defined in appendix The resulting expression for the propagator now has 
the form 

M M M 11 

S = So + So— So + So— So— So + ... = {So^ + M)-i (57) 

where M = M{p) is the momentum dependent mass defined above. There is just one 
thing to do: We have to solve the circular dependence 

& s 

nR -^m(^ 1 (H) 



but /i is just a number, which makes the solution very simple. Inserting (|56|) and (|57| ) 
into (|52|) one gets the following equation for /i 

^ = + / ^ p4lm + M,)^ ^^' + "^^"^ + ^-^^ ^^^^ 
which may be solved numerically for different current masses m. 
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Some Phenomenological Results 

In the chiral limit (|53) reduces to 



^^' = nnjd^ = annp' + 0{nl) (59) 

/i^ is proportional to n^i and thus M is proportional to ^/ur in contrast to a linear 
dependence on obtained from a naive density expansion. 

In the last expression the denominator has been expanded in the density and 

a = p-2 y (r^p2^'4(p) = 6.6 (60) 

is a universal number. For the standard values of and p one gets 

= G.Gurp'^ = (lOOMeV)^ , 
M{p = 0) = 7.7p^/^ = SOOMeV . (61) 



The exact solution of ( [59| ) which has been obtained numerically by iteration, differs from 
the leading density value by 15%: 

M(0) = 345MeV . (62) 

The momentum dependence of the quark mass is shown in figure |^. The mass m + M{p) 
may be interpreted as the mass of a constituent quark. At high energies it tends to the 
current mass, at low momentum chiral symmety breaking occurs and the quark gets its 
constituent mass M(0). Note that this is not a pole mass but a virtual mass at zero 
momentum squared. 



Let us now take into account a small current mass m formally of the order y^rTfj. The 
selfconsistency equation now reads 



l = - + ^ + 0{nn) . (63) 
P P 



Solving it for p leads to 



Po < P = ^m + ^^m"^ + pI < m + po (64) 

For the strange quark p is increased by a factor of 2: 

p{m, = 150MeV) = 200MeV (65) 

It is interesting that + M(0) remains to be SOOMeV. For zero momentum the increase 
of the current mass is just compensated by an equal decrease of the dynamical mass M(0). 

From the propagator one can obtain the quark condensate 

{^^) := lim tTcD{S{x) - So{x)) = nJ d^tioiSip) - So{p)) . (66) 
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In leading order in the density one gets 

= ^ = (G;:,Gn/32vrV ■ (67) 

This leads to the following condensates for u,d and s quarks: 

i{uu) = i{dd) = (250MeV)^ , (ss) = 0.5{uu) . (68) 
For heavy quarks there exists a similar relation 

^(^^) = (G;.Gr)/487r2m + 0{m-') , (69) 

which leads within 10% to the same value for the strange quark condensate. This nicely 
confirms the hypothesis that the strange quark can be treated as a light quark as well as 
a heavy quark. This hypothesis is used in heavy to light quark matching formulas. 

Large Nc expansion 

Consider now the case of Nf light quark flavors u,d, s, . . .. The discussion of the one flavor 
case in the previous sections can be copied up to the pairing and cancelation of graphs 
which contain more than one M/ (^). This is still true in the multiple flavor case but now 
both quark lines in (|55D must have the same flavor because Mj always connects quarks 
of the same flavor. So we have the theorem: "every Mj occurs only once for each flavor". 
From this point on the discussion of the one flavor case breaks down because there are 
now graphs contributing to the propagator containing quark loops. The simplest new 
contribution has the form 




Is this contribution small in some sense ? Yes it is ! Quark loops are suppressed by a 
factor l/Nc- In perturbative context this is extensively discussed in p8[, in instanton 



physics it was first used by Although 1/3 is not a very small number the large Nf. 
expansion seems to be a good approximation in various cases. 

Consider a graph and add to it a new quark loop consisting of 

N new instantons S > N 

S instanton scatterings ~ 

This multiplies the graph (see table 0) by a factor loop = 0{N^~^'^~^). The following 
cases are possible: 
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Parameters 


ur, p, Nc, (gR) 


Instanton density 
Instanton radius 
Number of colors 
Coupling constant 


n = N/Vi = urN, ^ (200MeV)^ 
p ^ (eOOMeV)-^ 
iVc = 3 

9 = Or/Nc (not used in semicl. limit) 


Gluon condensate 
Quark condensate 
Constituent mass 
Quark mass 
Gluon mass 

Ghost mass 
Meson correlator 


(GG) = {G^^Gnm^'' =■■ riRN, 
(ijij) ^ 0.39Ncp-^y/Wi ^ (253MeV)3 
M{p) ~ py/nE 

Mq^arkiO) ~ 7 .7 pn]{^ ^ 300MeV 
MgiuoniO) ~ 10.9 pn^/' ^ 420MeV 24 
MghostiO) ~ 7.7 pn]l^ ^ 300MeV|24| 


Quark loop 
Instanton scattering 
Instanton occurance 
Gluon loop 
Instanton scattering 


N-'pur'/' 

-<I) {N^ - 1)-' 



Table 1: Dependence of various quantities on the parameters of the instanton liquid 
model nR,p,Nc, (gR). 



1 + N~S 

= 1 <^=^ M = N <^=^ all instantons are new 

the loop is disconnected 
= <^=^ one old instanton <^=^ the loop is a tadpole 
< the loop is supressed by at least one l/Nc 



Disconnected graphs are cancelled by the denominator and tadpoles have been absorbed 
in ur. So quark loops are indeed suppressed in the limit iY^ — > cx). The same is true for 
gluon loops. This can be seen by the same argument using the Nc dependences from table 



4 Four Point Functions 



Introduction 

In the last section we have derived comprehensive Feynman rules within the zero mode 
aprroximation: 

' = ^ i"^- = ^o{p) 

P + tm 
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A quark of a given flavor can scatter only once at a instanton / via Mj. Tadpole graphs 
are not allowed, they are absorbed in the renormalized instanton density which has to 
be used when averaging over the instantons. In the large Nc limit dynamical quark loops 
are supressed and /i can be determined by (p8|). In the chiral limit 



/i^ = Q.QriRp^ + 0{n 



R) 



(70) 



In this section we want to calculate 4 point functions in the case of two quark flavors of 
equal mass in the limit ^ oo within the zero mode approximation: 



r' 



d{p-s + r-q)Urr'{p,s,q,r)= r 
dxdydzdw (.i{vy-iM-sx) (q | Tt/; (x) Ft/^ (y) F V (w) 1 0) 



(71) 



The il) fields are u ox d quark fields arbitrarily mixed. Without restriction to generality 
we have taken the correlator to be a color singlet . These 4 point functions can be used 
to study meson correlators (see chapter ^ 



n(t) = ^^KJ^^ 



(dpds) J {dqdr)Il{p, s, q,r) 
jr{x) = tjjTtp{x) 



dxe''^{0\Tjrix)jr'mO) 
jr'(0)=^rV(0) , 



(72) 



j{dpds)= jdp^sHp-s-t) , j{dqdr) = j^q^vHq-r-t) , 



t = p 



q — r 



Large Nc approximation 

The most general graph for the quark propagator is a sequence of different instantons M/, 
according to the rules above. Similarly the most general graph for the triplet correlator 
{{ud){du)) consists of two quark propagators each containing every insanton only once. 
However, the u and d propagator may contain common instantons, e.g. 

p q 




It is always assumed that the left and right hand sides of the graphs form color singlets. 
Non-common instantons can be averaged like in the propagator case to yield graphs of 
the form ^ ^ 

where the thick line represents the full propagator 
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^ = ( ' + ^^5^^ + + 

\ / /^Jt^., 

It can be shown that non planar diagrams Uke 



are suppressed by l/Ai'c, where again the left and right hand sides of the graphs have to 
form color singlets. So only the ladder diagrams shown in ( |75D contribute to the triplet 
correlator. 

One might think that the mixed correlator {{uu){dd)) is zero because the graphs neces- 
sarily include quark loops, or that only the two loop graphs contribute — but this is not 
the case! Let us first state the result and then discuss it. Graphs contributing to the 
mixed correlator are chains of quark bubbles 




(74) 

Application of the A*";, counting rules shows that this chain is of order l/N^. Taking the 
color trace at the left and right hand side of the chain we see that the mixed correlator is 
of order N^- Using it is clear that the triplet correlator is of the same order. 



What is wrong with the derivation of quark loop suppression in the last chapter ? The 
main assumption was that every graph containing a loop can be constructed from a graph 
not possessing this loop by simply adding the loop. Eliminating a loop from the bubble 
chain (|74D yields a disconnected graph, but we only consider connected 4 point functions. 
So the quark loop chain cannot be constructed in a way needed to prove quark loop 
suppression. 

In the case of meson correlators one can take another point of view. The disconnected 
two loop contribution is of order A*"^ but (except for the scalar case) the contribution is 
zero. So the bubble chain is a subleading graph of order N^. and nothing has been said 
about the form of subleading graphs. 

Nevertheless, all connected graphs can be obtained starting from (|7^) by adding further 
instantons and bubbles — but now Nc counting rules tell us that every attempt results in 
a 1/Ac suppression. Therefore the bubble chain is the most general leading order graph. 

Nothing has to be changed for the correlator {{dd){dd)) except that chain (|7^ must start 
with d. 
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To calculate the connected 4 point functions one must now average and sum up the chains. 
Alternatively this can be represented in recursive from usually called Bethe Salpeter 
equations: 



u 



G 



u 



d 



u 



Ih^ 

d 



+ 



u 



u 



d 



d 



G 



u 



d 



u 



u 



H 



d 



d 



u 



u 



d 



9 

d 



u 



u 



d 



d 



K 



d 



d 



(75) 



d 



K 



d 



d 




Solution of the Bethe Salpeter Equations 

Before solving the BS equations we have to construct the kernel. The l.h.s. of the kernel 
always forms a color singlet because of the restriction to color singlet correlators. 

Contracting the color and Dirac indices on the l.h.s. and using the formulas of appendix 
one gets 

r : ) = ^r/Mr)i^i{s)^^l^Mp)i^i{Q)(/)i= (76) 



s —h 



_f, — J. 



I {ifiy 



= -^p(f'{p)q(f'{q)r(f'{r)s(f'{s)d{p - s + r - q) {tiD [T 
The kernel can now be determined to be 



1 ± 75 \ 1 ± 75 



p-f- 



S —h 



-f-q 



-f—r 



V—h 



s —h 



-f—q 



-4 — r 



P 



nRN^ 



M,M,MrMJ{p -s + r- g)(5t^X + 75^75^)^^ 



(77) 



The result is just proportional to the nonlocal version of the 't Hooft vertex between color 
singlet states. 
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The solutions of the BS equations have a very similar structure: 



p-f- 




—f-Q 




A 




s —f— 




—f—r 



MpMqMrMsdij) - s + r-q) 



{7i 



Aq and A^ are scalar functions depending only on t=p — s = q — r. The proof is simple: 
The kernel has the structure ([78|) with Aq = A^ = 1. The product of two vertices yields 
the sames structure: 



—f— 


A 




B 


—f— 


—f— 


< 
— 


—/— 









s-f- 



AFF 



-f-q 



nRNc 



MpMgMrMsd{p -s + r-q) 



{AoFoBo{t)6X6l + A,F,B,{tH\jlrWi5X 



(79) 



Fo(t) = - ({dpds)—MpM,iiD{S{p)S{s)) , 
J riR 

F,{t) = - ({dpds)—MpMsiiD{S{p)i^S{s)^,) 
J riR 



^0) 



In other words, the vertices of structure (ffq ) build a closed algebra. The reason for 
this simple result is that the kernel is a simple product function up to the momentum 
conserving 5. 

Using (frzl), ([78|) and (^) the BS equations (|75|) reduce to primitive algebraic equations 
for Go/5(^), HQ/^{t) and i^o/sW: 

Ga/bif) = 1 + -^0/5(^)^0/5 (^) , 
Ho/5{t) = -1- Fo/5{t)Ko/5it) , 

KoM = -Fo/5(t)i/o/5(t) , 

with the solution 



G 



1-F 



H 



1 -F2 



K 



(81) 



where we have suppressed the index 0/5 and the argument t. 
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Triplet and Singlet Correlators 

Because of isospin symmetry SU{2)f, mesons form triplets and singlets. Replacing t/ji/j in 
(^) by the triplet and singlet combinations (borrowing the notation from the pseudoscalar 
correlator) 



TT 



one gets 



--j={uu — dd), 7i~^ = ud, TT = du, rj = --j=(uu + dd) , 



(82) 



^ 2 



u 




u 




K 




u 




11 



u 



u 



H 



d 







u 




H 


+ 


~1r 




u 



(j I 1 (j, 



d 



K 



d 



Therefore = K — H = jzp- This coincides with G = for the charged triplet 
correlator ((7r^)(7r^)) as it should be. In the singlet case we get C = K + H = —j^- 

When adding propagators in ( [78l) to the external legs the final result for the connected 4 
point function is 



U'^"^r{p,s,q,r) = '-JMpM.MrM, [Co{t)tTD{S{s)rS{p))tTD{S{q)rS{r)) + 



1 



Q>{t)tiD{S{s)TS{p)^,)tiD{S{q)TS{r)^,)]{M) 
+ for the singlet correlator 

(85) 



Fo/5(t) ±1 — for the triplet correlator 

Fo/^it) are defined in (pO|). The correlators of a singlet with a triplet current are zero as 
expected. 

The following graphs may contribute to the disconnected part: 



P 



r 



r = N^tiD{TS{p)T'{s))i{p - q)^{r - s) 



(86) 



P — I I — 9 



Y r' 

- r 



-N^tiD{rS{p))tiD{r'{q))^{p - s)5(g - r) 



depending on the flavor structure of the correlator. Note that the second two loop term is 
of the order N^. However, as discussed above, in most applications it drops out or yields 
an uninteresting constant or only the connected part is considered anyway. 

For the triplet and singlet case one gets 

for triplet 
2 • (86) for singlet 



d{p-s + q~r)U'^'^r{psqr) = NctTD{TS{p)r'S{s))d{p-q)d{r-s) + ' 



(87) 

The 4 point functions obtained in this section will be discussed in the following chapters. 
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5 Correlators of Light Mesons 



Analytical Expressions 

In the last chapter we have calculated various quark 4 point functions ([71|). The meson 
correlators or polarisation functions are just local versions of these vertices and can be 
obtained by simply setting x = y and z = w. In momentum space the meson correlators 
have the form 



Urr{t) = W^'^'it) + n'=°""(t) = ^K'y^= (88) 
'{dpds) lidqdr)n{p,s,q,r) = - ldxe''^{0\Tjrix)jr'mO) , 

Mxy/' = ^iuTu{x) - dTd{x)) , jriOr^' = ^iuT'u{0)-dT'diO)) , 

J {dpds) = J dpds d{p — s — t) , J {dqdr) = J dqdr d{q — r — t) , 

t=p — s = q — r . 

From the explicit expressions of the 4 point functions obtained in the last chapter one can 
get, up to integration, analytical expressions for the meson correlators. The following list 
is a complete summary of all formulas needed to evaluate the meson correlators: 

Uf^?{t) = N,J{dpds)tin{rS{p)T'Sis)) , 

n-T(t) = -N,{Co{t)rUt)TUt) + C5{t)rUt)TUt)) , (89) 

^ 1 + for singlet correlator 

Fof^it) ±1 — for triplet correlator ' 

r°(t) = -L [{dpds) JMpMstTDiSip)TS{s)) , 
TUt) = -L [{dpds) JM,M,tTD{S{p)TS{s)j: 
Fo{t) = — [{dpds)MpM,,tTD{S{p)S{s)) , 

F5{t) = — [{dpds) MpM,tTD{S{p)-f5S{s)-f5 



5 J 



pip'{p) = 27rpz—[Io{z)Ko{z) - Ii{z)Ki{z)]^=pp/2 , 



d_ 

dz 

2ip'^ {p)Mj 
p2 + (m + Mp)2 
urN, = (200MeV)^ , p = (600MeV)"^ 



p = m 



J + m + M„ r 
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Analytical Results 

Performing the Dirac traces leads to the following expressions: 

4 f MpMs{±{ps) - MpMs 



Fo/5{t) = fidpds 



4 r JMpMs{±{ps)-MpM, 
— lid'ods) — 



\ / If' R ^ 



4^ r jMpMs{MpSf,-MsPf,) 



/^J' ' {p^ + M^){s^ + M^) 
Mp = m + Mp 

All other vertices F are zero. Consider the one instanton vertex ( [77| ) (the kernel). It 
contributes only to the scalar and pseudoscalar correlator. From this observation one 
may have predicted that the connected part of all other channels is small because a 
contribution has to be a multi-instanton effect. Indeed, they are all zero as seen above 
except for the axial correlator. Due to an extra factor M ~ ^/nR in the numerator of F^g 
the connected part of the axial correlator is suppressed by 0{nR) therefore it is small as 
expected and will be neglected in the following. 

Furthermore we will restrict ourself to the chiral limit, taking m = 0. Using the selfcon- 
sistency equation (|90|) one can see that 



F,{t = 0) = l (92) 

is leading to a pole at t = in the pseudoscalar triplet correlator due to the -p5(t) — 1 
denominator in (|89|). This is the massless Goldstone pion one expects in the chiral limit. A 



more extensive discussion can be found in [|T9[. On the other side in the singlet correlator 



the minus sign is replaced by a positive sign and there is no Goldstone boson in this 
case. Thus, the (two flavor) rj' meson is massive ! Unfortunately we cannot make any 
reliable prediction of the rj' mass because the kernal is very repulsive in this channel and 
no boundstate is formed. There have to be other attractive forces, e.g. confinement forces, 
to built an rj' boundstate. Similar things happen in the scalar triplet channel (compare 
figure and |^). But the most important thing is that there is no massless pseudoscalar 
singlet meson which is an important step towards discussing the f/(l)A-pi'oblem. 

It is interesting to see that in leading order in the instanton density -fo(^) = ^-^5(^)5 which 
leads to a massless pole in the scalar singlet correlator. Numerically the cx-meson indeed 
turns out to be very light. The experimental situation is rather unclear. 



Spectral Representation 

To extract phenomenological information from the meson correlators we make use of the 
spectral representation 

//•oo 
dxe'P'-{0\rjr{x)jr'mO) = da^ D{a,x)p{a') (93) 
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where 

p(/) = (2vr)^5:5(p-g„)(0|jr(0)|n)(n|jf(0)|0) (94) 

n 

is the spectral density and 

D{m,x) = / ct^p^^, ^ = , „ ^ {mx)Ki{mx) (95) 

is the free propagator of mass m in coordinate representation. We have chosen the coordi- 
nate representation of 11 to be able to compare the plots directly with lattice calculations 
and with numerical studies of the instanton liquid 

The spectrum consists of mesonic resonances and the continuum contributions. If one is 
only interested in the properties of the first resonance one might approximate the rest of 
the spectrum by the perturbatively calculated continuum. 

One might think that the disconnected part only contributes to the continuum and the 
connected part will yield the boundstates. But this is not the case. On one hand, Bethe 
Salpeter equations have bound as well as continuum solutions. On the other hand consider 
a theory with weak attraction between particles of mass m. It is clear that there is only a 
cut above 2m and no boundstate pole in the free loop. However in the exact polarization 
function only a small portion of the continuum will be used to form a pole just below the 
threshold because the attraction is only weak. The Euclidian correlator will hardly be 
changed. Therefore, assuming weak attraction, we can already estimate the boundstate 
mass from the disconnected part. Of course in this example we need not calculate anything 
because we know that the bounstate mass is approximately 2m with errors of the order 
of the strength of the interaction. 

Assuming that all other forces neglected in QCD so far, especially perturbative corrections, 
are small and attractive in the vector and axialvector channel, we can obtain boundstate 
masses although in these channels up to our approximation there is no connected part. 
But things are less trivial than in the example above because the quarks do not posses a 
definite mass and we have to inspect the correlator to extract the meson masses. 

Let us start with the scalar and pseudoscalar correlator. The lowest resonance of mass 
is coupled to the current with strength 

A* = (0|ji/5(0)b) . (96) 

The rest of the spectrum is approximated by the continuum starting at the threshold E^. 
* means vr, t], 6 or cr (see table i?* is typically of the order 1.5 GeV and therefore 

the continuum can be calculated perturbatively. The spectrum thus has the form 

Pids) = ms - ml) + ^^e{s - E'^) . (97) 

Inserting p into (p3|) one gets 

nfj^x) = XlD{m,,x) + E{E,, x) , (98) 
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Correlator 






r = r' 


1=1 


1=0 


Pseudoscalar 


n5 = 


= (J5J5) 




TT 


7]' 


Scalar 


ni = 


= (iiji) 


Id 


6 




Vector 


n,M = 


~ {jfijfj) 


7/. 


P 


UJ 


Axialvector 






7m75 


ai 


h 



Table 2: Mesonic correlators 



E{E,, X) = ^ {2K^{E,x) + {E,x)K2{E,x)) ^ . (99) 

n^x" lb 7i^x° 

In the next section ^m^,, A* and E^, are obtained by fitting the phenomeno logical ansatz 

n-'^**(x) to the theoretical curve 11'^'"'™ (x) in the Euclidian region where the theoretical 

calculation is reliable. 

Consider now the vector and axial vector correlator. The vector current is conserved, thus 
the correlator is transverse and only the vector meson can contribute. In the chiral limit 
the same holds true for the axial current. In the singlet channel one has to be careful 
because there are two currents. A conserved one and a gauge invariant one which contains 
an anomaly. Up to now we have only calculated the correlator of the conserved current. 
Nevertheless to leading order in the instanton density the two correlators coincide and 
should be both conserved. 

For conserved vector and axial currents the spectral function is transverse: 

p^:^ip') = i-s,.+^)p'^\p') . (100) 

The coupling of the vector and axial meson to the current is given by 

^X.e, = {0\jf\0)\p) (101) 
where is the meson polarization. The spectral and polarization functions have the form 

= 3A^5(s - m^) + |^e(. - E,2) , (102) 
-Uj:^'^'\x) = 3XlD{m,,x) + 2E{x) . 
Here * means p, u, ai or /i (see table 0). 



Plot & Fit of Meson Correlators 

The meson correlators are shown in figure ^ - The numerical evaluation of the integrals 
is discussed in appendix p. The correlators are normalized to the free correlator 

n;/,M-i . K?--^ ■ (103) 

The diagrams therefore show the deviation from the perturbative behaviour. The meson 
parameters obtained by fitting the parameter ansatz to the theoretical curve are sum- 
marazied in table 0. Shuryak & Verbarshot plj have obtained the same parameters from 
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a numerical investigation of the instanton liquid model. Their values are also shown in 
table ^. The parameters of the vector channel coincide extremely well with ours. For the 
TT and a meson there is a large discrepancy in the masses but this is not surprising: We 
are working in the chiral limit thus the pion mass has to be zero. A similar argument 
holds for the a meson as discussed above. The couplings fit very well. The discrepancy in 
the axial channels can have various origins which are under investigation. Alternatively 
one may directly compare the graphs. They coincide very well even in cases where a spec- 
tral fit does not work very well like in the 6 and t]' channel. The conclusion is that the 
terms neglected in our analytical treatment, but included in the numerical study ||2T|, are 
small and usually give an correction less than 10%. These are contributions from nonzero 
modes and higher order corrections in l/Nc- This is again an example for the surprisingly 
high accuracy of the 1/Nc expansion. In the case of strange quarks the nonzero mode 
contributions will become more important. 

Finally, one should compare the numbers with experiment. As far as known, these num- 
bers are also listed in table 0. A general discussion of the meson correlators and compar- 



ision with experimental results can be found in |21] 



6 Conclusions 

Various QCD correlators have been calculated in the instanton liquid model in zeromode 
approximation and l/N^. expansion. The l/N^ expansion should be seen as a substitute 
for the density expansion which fails in the presence of light quarks. We have extended 
the work [|19] by including dynamical quark loops. In contrast to the original "pertur- 



bative" l/N^. expansion not all quark loops are suppressed. In the flavor singlet 
meson correlators a chain of quark loops survives the l/N^ limit, causing the absence of 
a Goldstone boson in the pseudoscalar singlet channel. The analytical expressions for 
the meson correlators have been evaluated numerically and meson masses and couplings 
have been obtained from a spectral fit. Except for the pseudoscalar singlet and the scalar 
triplet channel, the spectral ansatz matches very well the theoretical curve. The num- 
bers obtained in this way are consistent with those obtained by numerical studies of the 
instanton liquid within 10%. Comparision with experiment is also quite satisfactory. 
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A Zeromode Formulas 



The covariant derivative p in the background of an instanton has one zeromode 

(t^ - 4j)^Ijj = (104) 
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Meson 


JG(^JPC^ 


m^MeV] 


VK[MeV] 


E,[MeV] 


source 









508±1 


1276±33 


1/iVe 


TT 


i-(o-+) 


142±14 


510±20 


1360±100 


simulation 






138 


480 


— 


experiment 








? 


? 


1/iVc 


T)' 


o+(o-+) 




? 


? 


simulation 






960 


? 


— 


experiment 








? 


? 


1/iVc 


6 


i-(o++) 




? 


? 


simulation 






970 


? 




experiment 






433±3 


506±3 


1446±20 




a 


0+(0++) 


543 


500 


1160 


simulation 






? 


? 




experiment 






930±5 


408±4 


1455±33 


1/iVe 


P 




950±100 


390±20 


1500±100 


simulation 






780 


409±5 




experiment 






930±5 


408±4 


1455±33 






o-(i--) 


? 


? 


? 


simulation 






780 


390±5 




experiment 






1350±200 


370±30 


1050±80 


1/iVc 




1-(1++) 


1132±50 


305±20 


1100±50 


simulation 






1260 


400 




experiment 






1350±200 


370±30 


1050±80 


1/iVc 


/l 


0+(l++) 


1210±50 


293±20 


1200±50 


simulation 






1285 


? 




experiment 



Table 3: Meson mass m*, coupling constant A* and continuum threshold E^, obtained 
within the instanton liquid model in this work (1/-/Vc expansion), from numerical simula- 
tion and from experiment. 
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The solution for an instanton of topological charge Qj = ±1 in singular gauge in co- 
ordinate and momentum space is given in terms of modified Bessel functions Jq/i and 



Ko/i 



M^ + Zl) = y2<^(x)^X/ , fix) = „|^^|^^2^^p2)3/2 



TT 



^j{p) = ^f^^K 'P'^xf , vip)= Id^xe'^'^^ix) (105) 



y ■ PP ^ ^ 



2 ^ some color & Dirac column 
As usual only the projectors can be given in a covariant waj 



= (r, ±^) 



Xfx^ = ^(7.7.^)(f//T-;r±t/t) (106) 

xm = T'-{l,^-^mT^U\) (107) 

Uj/j e SU{Nc) is the orientation matrix of the instanton I/J. Although multi- instanton 
effects are studied in this work formulas like ( |10(j| ) concerning the overlap of different 
instantons are not needed. ( |106| ) is only needed for J = J in the following special cases: 

Nc{xfxf)i = tTDxfxf = ^(^^) (108) 

XX = tTcDXX = 1 (109) 
Taking space, color and Dirac part together we get 

tp\{p)Mp) = ^v'\p) =^ f(^4^\{p)Mp) = 1 (110) 



B Numerical Evaluation of Integrals 

The integral expressions for the meson correlators have been evaluated numerically. Two 
types of operations have to be performed: 

1. Convolution of Lorentz covariant functions (Fo/SjFr) 

2. Fourier transformation (FT) of the correlators to coordinate space 

Euclidian space the bar operation is the same as the hermitian conjugate: no 70 is introduced, 
thus X = X'- 
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Fourier Transformation 

Let us first consider tlie FT generalized to d dimensions: 



= F,{f,,...,J{x) = Jd^e-'^^f,,...M (111) 

For a scalar spherically symmetric function / = the FT reduces to a one dimensional 

integral 

poo / rn \ '^^'^ 

FMX\p\)}i^)= [Tm) fMJd/2-iHx\)\x\dm (112) 
Jo \27r|a;|y 

where are Bessel functions. 

If X is not too large and / decays rapidly, the integration can be performed with gaussian 
(or other) integration methods. If the decay is too slow one has to subtract the asymptotic 
part from / thus improving convergency. The FT of the asymptotic part can be performed 
analytically and has to be added to the numerical FT of the reduced function. 

The FT of a general Lorentz covariant function can also be reduced to (|112|) with (for- 
mally) an increased dimension d: 

Fd{pJ{\p\)}{x) = 2mx^F,+2{f}{x) (113) 
F,{p,pJ{\p\)}{x) = j^\{S,.-^)F,{p'f{\p\)}{x)- 



x^ 



(V - d^){AnF,^2{f}{x) - 47rVF,+4{/}(x)) 



Convolution 

Now we treat the convolution integrals 

f*9{p) = Jc^qf{q)-9{p~q) (114) 

This type of integral can be reduced to the FT discussed above: 

f*9ip) = F^'{F4f}-F,{g}} (115) 

This is a quick and easy method for evaluating convolution integrals. For the disconnected 
part of the correlators it has the advantage that in coordinate representation F^^ can be 
dropped. The disadvantage of this formula is that the FTs involve oscillating integrals 
which are numerically problematic. If the back-transformation F^^ is needed as in the 
case of Fo/5 and Fp it is better to perform the convolution directly. Similar to the FT the 
convolution can be reduced to the scalar case. The convolution of two scalar functions 
can further be reduced to a two dimensional integral: 

f*3{p) = 2j^,/ii^^!^2)\ L drldef{r)g{^p^-2\p\rcose + r^){rsmef-'r (116) 
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It is again evaluated with gaussian integration methods. The second advantage is that 
there are no problems with slowly decaying functions. Sometimes there are large cance- 
lations between different terms. In this case it is essential to use nonadaptive integration 
methods because they will not result in a loss of accuracy. 

The explicit reduction of the various correlators to the basic forms ( |112| ) and ( |116| ) is more 
or less trivial. The selfconsistency equation has been solved by iteration. The results are 
plotted in figure H - B 
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ifigl3.pic 



Figure 1: 5* is minimal at the river, small in the valley and large in the mountains, 
Therefore Z = Jdxidx2e~^^^'^'^^^ is dominated by the valley. 
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Figure 3: Pseudoscalar triplet correlator normalized to the free massless quark correlator. 
The pion coupling constant and the continuum threshold E^^ are fitted in order to 
match the spectral ansatz with the theoretical sum of the free and the conrzected part. 
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Figure 4: Pseudoscalar singlet correlator normalized to the free massless quark correlator. 
There is a strong repulsion in this channel and no boundstate is formed. The theoretical 
curve is compared to a curve obtained from a pure continuum spectrum above E^,' . 
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Figure 5: Scalar triplet correlator normalized to the free massless quark correlator. There 
is a strong repulsion in this channel and no boundstate is formed. The theoretical curve 
is compared to a curve obtained from a pure continuum spectrum above Es- 
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Figure 6: Scalar singlet correlator normalized to the free massless quark correlator. The 
a mass mg. and coupling A^. and the threshold E^j are obtained from a spectral fit. 
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Figure 7: Axial vector correlator normalized to the free massless quark correlator. The 
triplet and singlet correlator are equal because the connected part has been neglegted. 
The Oi and fi mass, coupling and threshold are obtained from a spectral fit. 
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Figure 8: Vector correlator normalized to the free massless quark correlator. The triplet 
and singlet correlator are equal because the connected part is zero. The p and lj mass, 
coupling and threshold are obtained from a spectral fit. 
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